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Abstract 

We revisit the recently proposed method for the determination of the vrvr scat- 
tering length combination ao — a2, based on the study of the spectrum in 
K+ — > 7r"'"7r''7r'^ . In view of a precision measurement, we discuss here the ef- 
fects due to smaller absorptive contributions to the 7r'''7r'^7r^ and 
7r^7r^7r~ amplitudes. We outline a method of analysis that can lead to a precision 
determination of ao ~ ^2 and is based on very general properties of the S matrix. 
The discussion of final-state rescattering and cusp effects is also extended to the 
two Kl — > Svr coupled channels. Thanks to the present work, the theoretical error 
on the ao — 02 combination extracted from the — > 7r"'"7r^7r*^ spectrum is reduced 
to about 5%. A further reduction requires the evaluation of the effects arising from 
radiative corrections. 



1 Introduction 



In it was shown how rescattering of the final state pions in Tr+Tr^vr" produces 

a prominent cusp^ in the total energy spectrum of the vr^vr^ pair, whose amplitude is 
proportional to the — a2 combination of the tttt S-wave scattering lengths. 

The combination oq — 02 is a very interesting quantity: a benchmark observable to 
determine the structure of the QCD vacuum and one of the few non-perturbative pa- 
rameters which can be predicted with excellent accuracy from first principles. Recent 
calculations [4 , that combine Chiral Perturbation Theory (CHPT) [3j with Roy equa- 
tions ISIEI, lead to the precise prediction (oq — 02)^^+ = 0.265 ± 0.004. So far, this high 
theoretical precision has not been matched by a similar experimental accuracy. The best 
direct information on tttt scattering lengths is the one extracted from Ke4, decays by the 
BNL-E865 experiment [Zj, which is affected by a sizable (~ 6%) statistical error. Given 
the intrinsic statistical limitation of Ke4 decays with respect to the dominant K ^ Sir 
modes, and the different nature of systematical errors (including the theoretical ones) 
involved the extraction of vrvr scattering lengths in the two cases, it is definitely worth to 
explore in more detail the proposal of Ref . jT] . 

The NA48 runs of 2003 and 2004 have produced ^ 10* 7r+7r°7r° decays, of 

which a few millions are in the tt'^tt' threshold region with excellent St^tt resolution. Since 
the cusp induced by the tttt rescattering is a ~ 10% effect on the 7r°7r° spectrum, at a pure 
statistical level it should be possible to determine its amplitude with ^ 1 — 2% accuracy. 
In order to extract from this measurement a value for oq — 02 with a similar precision, it 
is necessary to reduce the theoretical uncertainties of the simple analysis proposed in pp. 
The present paper is a first step in this direction. 

As already noted in P], the procedure presented there was incomplete for three main 
reasons: 

1. It did not take into account the effect of radiative corrections. 

2. It omitted higher terms in v, e.g. terms in the imaginary part of the amplitude. 

3. It omitted contributions from higher order rescattering effects. 
^ The existence of this cusp was first discussed in 
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The effect of radiative corrections to the 7r+7r°7r° decay in general, and to the 

amphtude of the cusp in the vr^vr^ spectrum, will not be discussed here. The radiative 
corrections to K —>■ in decays have been evaluated (see e.g. Ref. jH|) to be of the few 
percent level, and dominated by Coulomb corrections. We expect radiative corrections to 
the cusp amplitude to be not larger than this. 

The second omission is minor, as the value of ao — 02 is given by the term proportional 
to V, while the term in can be introduced as a free parameter in the experimental fit, 
and its possible prediction in CHPT is of lesser importance than that of the scattering 
lengths. The evaluation presented here includes these effects. 

In this paper we will concentrate on correcting the third omission. We will show 
how the unitarity and analyticity of the S matrix elements can lead to a systematic 
expansion of the —>■ tt+tt^tt^ and tt+tt+tt^ amplitudes in powers of the tttt 

scattering lengths. The usefulness of this expansion derives from the relative weakness 
of vrvr scattering that, in turn, is a general consequence of the pseudo-Goldstone-boson 
nature of the pions and of the smallness of light-quark masses (or, in one word, a general 
consequence of CHPT). Rescattering effects in K 3tt decays have already been widely 
discussed in the literature in the framework of CHPT (see e.g. Ref. [Hll^)- However, 
previous analyses have been performed only up to next-to-leading order in the chiral 
expansion and - with the exception of Ref. (TU] - ignoring isospin breaking effects. The 
approach presented in this paper is more general and particularly well suited to discuss the 
cusp amplitude: we shall use the effective field theory only for an explicit estimate of the 
irreducible Svr Svr rescattering (that turns out to be negligible). For completeness, we 
shall also present a general parameterization of rescattering effects and cusp amplitudes 
in Kl —>■ Stt decays. 

The paper is organized as follows: in section |21 we shall introduce the definition of tttt 
scattering lengths used in the rest of this work, and we shall recall some basic properties of 
the S matrix. Section IHl is devoted to analyse the consequences of unitarity and analyticity 
on the structure of various tttt tttt amplitudes. In section 0] we shall present the 
systematic expansion of K ^ 3tt amplitudes in powers of the vrvr scattering lengths up to 
0{af), we shall also briefly discuss possible strategies for the data analysis. The results 
are summarized in the conclusions. 
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2 TTTT Scattering 

2.1 Two pion states 

Consider the S matrix element: 

{c,Pc;d,Pd\S\a,Pa;b,Pb) = {c,Pc;d,Pd\a,Pa;b,Pb) + iS'^iPf - Pj) 



The normahzation of the states is chosen as 

(C, Pc, d, Pd\a, Pa, b, Pb) = Sea Sdb S^{Pc - Pa) S^{Pd " Pb) + ^cb ^da S^{Pc " Pb) S^{Pd ' 

This normahzation is compatible with the field theoretical definition: 

\a,Pa;b,Pb) = ai{pa)al{pb)\^) 

[aa{Pa): OblPft)] = SabS^{Pb-Pa) 

If we change variables to total and relative four-momentum, 

P = Pa+Pb; k= {pa - Pb)/2 P' =Pc+ Pd, k' = {pc - Pd)/'2 
S'iPc - Pa) S\pd - Pb) = S\P' - P) 5%k' - k) 
SHpc - Pb) S'ipd - Pa) = S'{P' - P) 5'{k' + k) 

— * 

we can define the S-wave state in the center of mass, P — 

\P, q, a, b) = — ^ / d^k S{q - \k\) \P, k, a, b) 
Vinq J 

and verify that the normalization is 

{P', q', c, d\P, q, a, b) = S\P' - P)S{q' - q) {ScaSdb + ScbSda) 

2.2 Isospin states 

We will adopt a phase convention, inspired by a field theoretical treatment, where 

r|7r+) = -V2|7r°); r|7r°) = v^lO 



We note that this convention is different from that used in early I-spin analysis of K ^ Stt 
decays, e.g. in ref. jTTj and ^12J, but coincides with the one adopted in CHPT studies of 
these decays, e.g. ref. ^H]- For Is = we then find the three states 

|2,0) = '—^ ^ ^ (11) 



11.0) = (12) 
10.0) = + I" + (13) 



And for Jc. = 1, 



12.1) ^ ^II^^^ (14) 

II. i> ^ "^-°>;,"°-^> (15) 

These states are normalized as (see eq. El but note that the 1 = 1 states vanish for 
S-waves) 

(P', q', /', I'^\P, q, I, I,) = 25rj 6j>j^ 6%P' - P) 6{q' - q) (16) 

2.3 Low energy scattering and scattering lengths 

S-wave scattering means that the /j defined in ((H) does not depend on the direction of 
the relative momentum /c, but at most is a function of the CM energy E or momentum 
q. We than find easily that, working in the CM. frame (P = 0), 

(P',g',/|(S - l)|P,g,^) = A7:iqq'5{E' - E)5\P' - P)-^#= (17) 

-^-'(liT|)*<''--)^''^'-^'Jfe (^«) 



TT I 



iE[ + E 



5{q' -q*)5\P' -P)Mf, (19) 



with q* the CM. momentum required by energy conservation. In the non relativistic 
limit, 

(P', q\ /|(S - 1)|P, q,^) = Amqf,5{q' - q*)5\P' - P)^= (20) 

= n^6iq'-q*)SHP'-P)Mf, (21) 
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Neglecting 7t~^7i^ mass differences, fi = 171.^+ /2. For exact I-spin (that only makes sense in 
the limit of equal masses) we must have near threshold: 

S\P,q,I,h)=exp{2iqaj)\P,q,I,h) (22) 

so that 

{P', q', I', J^|(S - 1)|P, q, I, /s) ^ Atqaj 6ri 5i^i, 5%P' - P) 5{q' - q) + 0{q^) (23) 
Comparing now ()2()|) and ()23j) we find that near threshold 

Mu^^-^ (24) 

TT 

Identifying m^r with we will thus define: 

00 00 -D \A 8aoo"^7r+ / + _ \A\ I-spin 00 + 202 . s 

TT TT 7T Tc ReAloo = [TT^TT threshold) aoo — ^ (25) 

TT 3 
+ +0 \A 8a+o'^7r+ / + U UN I-spin ^2 

TT^TT — > TT^TT ReAl+o = [71^71 thrcshold) a+0 — ^ 17 (26) 

TT 2 

+ - T> \A 8axm^+ / + - -i 111 i-spin cto — ^2 , , 
TT TT n 71 tie Aix = [71 71 threshold — > (27) 

77 3 

+ - + - 15 8a+_m^+ / + - 1 1 IN I-spin 2ao + (32 /^.^x 
TT^TT ^ TT^TT KeAl+_ = (tt^tt threshold) a_|__ — > (28) 

TT 6 

7r+7r+ ^ 7r+7r+ ReA<++ = ^^±±^ (7r+7r+ threshold) a++'^a2 (29) 

TT 

For each process we have noted the threshold at which the scattering length is defined, 
and the value it would have in the limit of exact I-spin symmetry. 

The problem that must at some time be faced in comparing the result of cusp studies 
to the CHPT prediction for — 02 is that of taking into account radiative corrections. 
Note that the threshold region is one where I-spin is maximally broken. We will take the 
point of view that the quantity introduced in eq. ()27|) should be taken as a definition 
of the effective scattering- length combination measured from the cusp effect. The experi- 
mentally determined value for this quantity should be compared with a CHPT prediction 
which includes the effects of radiative corrections and I-spin breaking due to 7^ m^. 
The evaluation of these subleading effects can be subdivided into two separate tasks: 
computing their impact on the CHPT predictions of the various vrvr tttt amplitudes in 
eqs. and determining how they would affect the decomposition of the K ^ 37i 

amplitude presented in this paper. The first point has already been partially addressed 
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in the literature [T3] - it turns out to be only a few percent correction in the case of 
(^x ^Sj ^ and can be easily implemented in our decompostion. However, at the moment 
we are lacking of a consistent description of the second point, or the evaluation of I-spin 
breaking effects in the relation between tttt and K ^ Sn amplitudes. 

In the following we will proceed using eqs. (PH|) - ()29|1 as a definition of the different 
scattering length combinations. We shall use their expressions in terms of ao, 02 only as 
a first approximation, pending a consistent evaluation of all the I-spin breaking effects. 
Note the use of in these definitions, and of the vr+vr" threshold^ in (j2Sl)- 

In the case where the scattering occurs well above threshold, eqs. (I25|) - (j29|) should be 
modified to take into account the non trivial kinematical dependence of Aifi- Expanding 
up to linear terms in the kinematical variables s, t and u, we can neglect all higher modes 
but the P wave. The generic matrix element takes the form^ 

= !!!^zi^ + (30) 

TT 71 + 

with 

[■5 "^threshold) 



The ttij are the combinations of constant S-wave scattering lengths defined in eqs. fl25j) - 
((2ni), while the define the corresponding effective ranges. In the isospin limit, we 
can express the Vij in terms of the effective ranges of oq and 02, following the isospin 
decomposition reported in eqs. (f^ - (pn|) . According to the detailed analysis of Ref. 
these are given by tq = 1.25 ± 0.04 and r2 = 1.81 ± 0.05, values which are consistent 
with those recently reported in Ref. and also not too far from the lowest-order CHPT 
predictions r^^^ = 8/7 and r^^ = 2. 

The only two channels with non vanishing afj are the tt+tt^ tt+tt^ and tt+tt*^ tt+tt'^ 
ones. In the I-spin limit 

a+_ = a+o = ai/2 , (32) 

while the lowest-order CHPT prediction is af'^ = m^/(127r/^). 

^ At this threshold the cusp correction to the tt^tt^ — + tt^tt*' scattering amphtude vanishes. 
^ This expression does not include the effects of threshold singularities, whose structure will be dis- 
cussed in section 13 
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2.4 Cluster decomposition and the Operator notation. 

The S-matrix elements can in general be cluster-decomposed*^ into the sum of a connected 
part (in perturbation theory this is the sum of the connected diagrams), and one or more 
terms that are the product of connected terms, and correspond to the separate interaction 
of non overlapping subsets of the initial particles to yield non overlapping subsets of the 
final particles. Among the disconnected terms there may be some where one or more of 
the initial particles propagate without interacting at all. 

It will be convenient to express the S and T operators in terms of creation and 
annihilation operators for asymptotic states, so that we can write T as a sum of operators: 

S = 1 + ^T; T = ^T^,„ (33) 

Each of these operators can be expressed as: 

T^,. = E / n [d'p] ^\Pf - p-^< ■ ■ ■ • • • ^^-^m, ^^^^ 

where the sum is over particle types and the integral is over the three-momenta of the 
initial and final particles. We note that each Tm.„ can contribute to an n — > m transition, 
but also toanra + A; — > m + k transition, with k particles passing through without 
interacting with the others. In general the Tm,n term will contain a connected part, T^^, 
where the corresponding T/j contains a single S'^{Pf — Pi) factor, and other terms with 
two or more such momentum conservation factors, and can be expressed as well ordered 
products (annihilation operators on the right) of "smaller" T^, For instance in the 
case m = n = A we would find 

iT4,4 = ^Tj4+ : {tT2,2f : (35) 

In the case we will be interested in, the — > tt+tt^tt*^ and tt+tt+tt^ decays, we 

will be working with T^, „ terms that coincide with their connected parts. We do not 
thus need to explore the disconnected parts of Tm,n in more detail here. 
* For a discussion of cluster decomposition, see e.g. ^|^], and |19l Vol I, Ch. 3]. 
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2.5 Time reversal symmetry 

We will neglect the effects of time reversal and CP violation on — > Stt decays. We have 
very strong experimental limits on these effects, and the theoretical expectation is even 
smaller. Time reversal symmetry implies the relation 

{B\S\A) = {At\S\Bt) (36) 

where \Bt), \At) are the "time reversed states", that for pseudoscalar mesons amounts 
to changing the sign of all momenta, p —p. In the case oi K 3tt and tttt ^tttt , 
that arise in the following discussion, we can change the sign of all momenta with a 
combination of Lorentz transformations and rotations, so that we have simply 

{B\S\A) = {A\S\B) (37) 

Because of parity conservation, the same condition holds for the Sir — > Sir strong re- 
scattering. So that, neglecting CP breaking effects, S is symmetric for all the cases of 
interest for this analysis. 



3 Unitarity, analyticity and the tt+tt threshold 

Due to the presence of the square-root singularity connected to the 7r"'"7r~ threshold within 
the phase space for — > tt+tt^tt^ , we have to distinguish the two zones above and below 
the TT+TT" threshold. We can write the amplitude M. fi above the threshold in the form: 



Mfi^A + B\l- ^ S3>4m2+ (38) 

^3 

where both A and B are regular at the tt+tt" threshold. This expression can be analytically 
continued below the threshold, where it becomes 



Mfi^A + iB\l^^ ss<Aml+ (39) 

S3 

Applying unitarity above the threshold we can determine the imaginary parts of both 
A and B. Also, unitarity below the threshold determines the real part of B. The experi- 
mental data can then be analysed with the following procedure: 
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1. Parametrize the real part of A as a polynomial in the three independent kinematical 
variables, si, S2, S3, as outlined in (20] • 

2. Parametrize the tttt scattering amplitudes in terms of the scattering lengths and 
possibly additional parameters. 

3. Use unitarity to derive B and the imaginary part of A. 

It is best to work out the consequences of unitary in the operator formalism, where 
we can express the S operator in terms of the hermitian and anti-hermitian parts of the 
T operator, 

S = 1 + 2(R + il) then unitarity implies (40) 

21 = + or, solving for I, (41) 



1=1 — Vl — and as a power series in R, (42) 

I = -R2 + 1r^ + -R6 + —R\.. (43) 
2 8 16 128 ^ ^ 

Time reversal invariance implies that S is symmetric, so that the matrix elements of R and 
I correspond directly to the real and imaginary parts of T matrix elements. 

The last equation allows for a systematic computation of the imaginary parts in terms 
of the real parts of the scattering amplitudes. The utility of this expansion derives from 
the assumed smallness of the tttt scattering lengths. In the case of —>■ tt^tt^tt^ the 
first term in the development yields terms oc and oc af, and higher, while the further 
terms, starting with R^, will contribute corrections oc af and higher. 

3.1 7r*^7r*^ scattering 

Let us apply the ideas outlined above to tt'^tt^ scattering. We will work in the threshold 
region, so that we can neglect higher partial waves and any dependence of the amplitude 
M.00 on the t,u variables. We will also neglect higher (e.g. in) cuts, and only retain the 
first term in eq. so that we will neglect terms 0{af) and higher. We will also use 

as unit of energy, so that we write e.g. aoo instead of aQom^+. 
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Let us start by defining the "velocities" , 



vAs) = \l ^' (44) 



Ms) = ^^-^^^ (45) 
We can then write, in analogy to eqs. (jHHj) . 

A^oo = ^0 + Boov±{s) s > Aml+ (46) 

A^oo = ^0 + iBoov±{s) s < 4m^+, (47) 

and, for tt+tt^ ^ ti^tt^, 

M^ = A^ + B^v±{s) s > Aml+ (48) 

where Aqq, Bqq, A^, are regular at the tt'^tt" threshold. We can express Re (A) as a 
polynomial in s. We can simply write: 

Re (Aoo) = ; where aoo(4m^+) = aoo (49) 

71 



and similarly for tt+tt vr^vr^. 



Re {A^) = ; where a,(4m2 + ) = a,. (50) 

TT 

The 7r"'"7r~ intermediate state contributes to Im A^oo only above the 7r~^n~ threshold, while 
the 7r°7r° state contributes both above and below, so that we find 

ImA^oo = jv^is) (ReM.f e{s - Aml+) + ^voo{s) {ReMoo? + O(R^) (51) 
4 o 

where O(R^) indicates the neglected higher terms in eq. (021) • Evaluating eq. (|3T|) above 
the TT+TT" threshold, and neglecting terms O(R^), this translates into 

Im Aoo = ^[(Re Aoo)^ + ^^^(ReSoo)^] + ^ReA^ReB^-^^ (52) 

8 s 2 s 

7T S — AlVT? _l 

Im^oo = -[(ReA,)2 + ^(ReS^)2 + ^;oo(s)Rev4ooRe5oo] (53) 

4 s 
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and evaluating it below the tt+tt threshold 

ImAoo = ^^[(Re Aoo)^ + ^^^^^(Imi?oo)1 (54) 

Re Boo = -^^^^Re Aoolm Boo (55) 

From eqs. . we see that lmi?oo = O(R^) and Re -Boo = O(R^), so that, comparing 
eqs. and (fH^ we conclude that also Rei?2,. is at least ©(R'^). Neglecting terms of 
O(R^), the final result is 

ImEoo = jiReA^y = - (a,(s))' (56) 

4 71 

T) R ^^"^00(5) , .2 32t;oo(s) , , ..2 

ReBoo = — — ReAoo(ReA^) = aoo(-s)(a^(s)) (57) 

ib 71 

ImAoo = — 3 — [ReAoo] = (aoo(s)) (58) 

Re 5^ =0(R2) (59) 

3.2 TT+TT^ and 7r+7r+ scattering 

In the following we will also need expressions for tt+tt^ and 7r+7r+ scattering. The situation 
here is simpler, since at O(R^) there is only one intermediate state. However, in this 
case we are interested in kinematical configurations where the amplitudes are not close to 
threshold and P-wave contributions cannot be completely neglected. Expressing the latter 
in terms of the 1=1 scattering length in units of 171^^+ (we can safely neglect I-breaking 
corrections in this case), the real part of the amphtudes can be parametrized as 

Re{M,o) = + (60) 

Re(-M++) = (61) 

TT 

We have adopted the standard notation t = {pi — p'lY and u = {pi — p'^i^ ■, where 
7r+(]9i)7r°(p2) — ^ ^^(K) + ^°(P2)- The imaginary parts can easily be computed, but 
are not needed in the following. 

3.3 TT+TT scattering 

In this section we will consider tt+tt^ tt+tt" scattering, that enters in rescattering cor- 
rections to the TT+TT+TT^ decay amplitude. We will again work close to the threshold 
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region, neglecting higher partial waves and the kinematical dependence from t and u vari- 
ables. Here we meet with a new problem: in the case of vr^vr'^ tt^tt^ scattering we were 
able to apply unitarity below the vr^vr" threshold, and this was used to derive a value for 
Rei?oo; eq. (j37j) . In the present case moving below the tt+tt^ threshold implies an 

analytic continuation to an unphysical region. We will proceed to do this by considering 
a continuation in the 7r+ and tt^ masses, a procedure that is certainly legitimate in a field 
theory, such as CHPT, where the 7r+ — 7r° mass difference can be changed by introducing 
an extra mass term in the Lagrangian. We will then work out the consequences of unitar- 
ity in a situation where m,rO > "^77+ , and analytically continue the results to the situation 
where the masses have their physical value. 

Assuming now that m^o > 771,^+ we must distinguish the case where s > 4m^o, where 
both TT+TT^ and tt'^tt'^ can appear as intermediate states, and that where s < 4m^o where 
only the tt+tt^ state can contribute. We can write 

M+- = C± + D±voois) s>Amlo, (62) 

M+- = C± + iD±voo{s) s < Amlo, (63) 

As in the case of the 7r°7r° scattering we start by defining the real part of C± in terms of 
the scattering length, 

71 

and applying unitarity at O(R^) we obtain: 

ImA^+_ = lvoo{s)(ReM.fQ{s-4mlo) + ^v±{s){ReM+-f + 0(R'') (65) 
o 4 

Evaluating eq. both above and below the 7r°7r° threshold, and neglecting terms 

O(R^), this translates^ into 

lmD± = l{ReA,f = -{a,{s)f (66) 
ReD^ = -^^^^ReAooiReA^Y = -^^^a+_(s)(a.(.))^ (67) 
ImC^ = ^{ReC^f = 1^ {a^4s)f (68) 

4 TT 



We omit the intermediate steps that follow the lines of the preceding section. 
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The continuation to the physical values of the 7r+ and vr" masses is simply achieved by 
using in the above expressions the correct masses for v±,voo and the physical values for 
the scattering lengths. 

4 K ^ Sn decays 

In the following we shall apply the results of the previous section to describe rescattering 
effects in — s> Stt decays. Our main interest will be on the ti^-k^ti^ channel, 

where the cusp effect is most prominent and useful for the determination of (oq — 0,2)5 
but we shall also consider the tt+tt+tt^ decay, whose amplitude is needed for the 

— s> 7r+7r°7r° analysis. Similarly, we shall discuss the Kl Svr^ decay, where smaller 
cusp effects — still proportional to (ao — 02) and related to the tt+tt^tt^ — > 3tt^ 

process — should also be visible. As in the previous section, we shall consider rescattering 
effects only up to 0(af) corrections to the leading amplitudes. To be more precise, we 
shall evaluate the full imaginary parts of the amphtudes at 0{ai) and the corresponding 
0{af) corrections to the real parts. 

Similarly to the vr'^vr'^ — 71^71^ case, we can decompose Tr+vr'^Tr'^ and — >■ 

TT+TT+TT^ amplitudes into a regular term and one that is singular at the vr'^vr" threshold. 
We will use the standard kinematical variables Sj = {px — Pwi)^, i = 1 ... 3, as specified 
in the Particle Data Group Review 1201, index "3" referring to the odd pion (vr^ or 
7T~ for the two decays). In particular, for K'^ -k^ti^tc^ , S3 coincides with the square of 
the CM energy of the tt^tt^ pair. We will thus write 

M00+ = A00+ + Boo+v±{s3) S3 > Aml+ , (69) 

M00+ = A00+ + iBoo+v±{s3) S3 < 4m^+ . (70) 

For the tt+tt+tt^ amplitude it will be convenient to separate the terms which 

contain the singularity in Si,S2 associated with the ti^ti^ threshold, and write: 

= + D^^l_voo{si) + Df_,^voo{s2) , (71) 

where Bose symmetry implies that 

■M(i){si, S2, S3) = A^(j)(s2, si, S3), (i = + H — , 00+), that translates into (72) 
-D+|_(si, S2, S3) = Z:'+|_(s2, Si, S3), etc. (73) 
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yloo+, -Boo+)C'++-) and -D++- are expected to be analytic functions of si, 82,83 in the 
physical region for the two decays, with square-root singularities at the borders associated 
with different vrvr thresholds. 

Unitarity will allow us to express Aioo+ aiid A4oo+ in terms of Rey4oo+, ReC++_ 
and TTTT scattering lengths. In the case of Rey4oo+, ReC++_ we adopt a parametrization 
inspired by the PDG tables, namely 

ReAoo+(si,S2,S3) =/?°(s3) = A° [l + (^0(s3-s°)/2m2+ + /iO(s3-s°)V2mt+] (74) 
ReC++_(si,S2,S3) =i?+(s3) = [l + (?+(s3-4)/2m2+ + /i+(s3-4)V2mt+] (75) 

where 

^ I 3s+ = ml, + 3m2, (K+ ^ tt+tt+tt- ) ^ 

The PDG tables also include terms proportional to (si — ^2)^, but their coefficients are 
small and compatible with zero^. They can be reintroduced, if needed for fitting a precise 
data-set, and one could also introduce higher powers of (53 — Sq). We can compare this 
parametrization with that adopted in the PDG if we neglect the other contributions to 
the decay amphtude discussed in this paper, namely the imaginary parts of y4oo+,C'++_ 
and the whole of i?oo+7-D^'^_, that give smaller contributions to the decay rates, except 
in the cusp region of tt+tt^tt^ . We then obtain: 

^7°'+ ^ ^7poo; /^°'+ ^ /ipoG - (^7°•+/2)^ (77) 

and, using the PDG average values, 

^° ^ 0.638 ±0.020; ^ -0.051 ± 0.01 (78) 

^+ ^ -0.2154 ±0.0065; /i+ ^ 0.0004 ± 0.004 (79) 

Interestingly, the PDG values suggest a quadratic term that is vanishing small in the 

—>■ TT+TT+TT^ amplitude. The small and negative value in the tt+tt^tt^ amplitude 

could simply arise from the effect on the previous fits of an undetected cusp in that decay. 

In this situation it would appear that the quadratic terms in eq. (fTIj) could be dropped 

^ The coefficient of the (si — 82)'^ term m |Aloo+P is known to be very small: it is listed in the PDG 
tables 20 as fc = 0.004 ± 0.007. The k coefficient for if + tt+tt+tt" is also compatible with zero, but 
with a slightly larger error. 
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b) c) 

Figure 1: Examples of K —>■ 3tt rescattering topologies at the two-loop level: a) single- 
channel nil scattering; b) irreducible Svr —>■ 3n contributions; c) 3n —>■ 3n amplitude due 
to multi-channel Tin scattering. 

at least in a first analysis that takes into account the cusp effect and other absorptive 
contributions. 

4.1 Three-pion scattering and two-pion cuts 

Our final goal is the evaluation of rescattering effects - and particularly the determination 
of the cusp amplitude - in the three-pion states produced by K decays. In general, in 
the case of Svr states, we can distinguish two basic contributions to the unitarity relations 
generated by ()43|) : those arising from rescattering of a pair of pions in a given channel 
- with a third spectator pion - (see e.g. fig. [T^) and those due to Svr 3it connected 
diagrams. At the level of approximation we are working, it is also convenient to distinguish 
between Svr Sir one-particle-irreducible diagrams (fig. [Dd) and Sn Svr reducible 
amplitudes due to multiple vrvr scattering in different channels (fig. [It). 
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The Svr Svr irreducible contribution is the only one that cannot be expressed in terms 
of vrvr scattering lengths, but it turns out to be safely negligible. A simple and reliable 
estimate of its size can be obtained using the lowest-order CHPT Lagrangian to evaluate 
the Stt —>■ Stt irreducible amplitude, and employing the non-relativistic approximation for 
the Svr states. In this limit, the irreducible scattering leads to a constant imaginary part 
of O{10~^). For instance in the —>■ tt^tt^tt^ case we find 

(ImAoo+)3. = -T^-f^(ReAoo+) ^ -4 x 10"^ (Re A00+) (80) 

where Q is the Q-value of the ti^ti^tt^ decay, = 130.7 MeV is the pion decay 

constant, and (Reyloo+) is the average of the real part of the amplitude over the Dahtz 
plot. This contribution, and others of a similar size also in the other channels, appears 
to be safely negligible at the 0(10^^) accuracy for the decay rates we are aiming for. 

The evaluation of the single-channel vrvr scattering in i^' ^ Svr decays proceeds exactly 
as for the Ain-K amplitudes discussed in the previous section. To this purpose, it is useful 
to observe that all the previous results can be recovered in a diagrammatic framework by 
considering the absorptive two-pion cuts of appropriate Feynman diagrams. In particular, 
the O(R^) contributions to the real part of the cusp amplitude (such as the expressions for 
Re -Boo and Re D± discussed before) can be derived by considering the two s-channel cuts 
of diagrams similar to the one in fig.[T^ (i.e. setting both the {vr^, tt^} pair and the {vr", vr^'} 
pair on shell). As we shall illustrate in more detail in the next section, this observation 
allows us to evaluate in a simpler way also the effects of vrvr scattering in different channels 
(fig. ^) and, in particular, to express them in terms of the vrvr scattering lengths. 

4.2 Rescattering in tt+tt^tt'^ 

We start by considering the rescattering of the two-pion pair leading to the final tt^tt" 
state. In the term of eq. ()43p the contribution of intermediate 7r"'"7r~ states — fig. |21 
— is given by 

(ImA^oo+)^+^- = / ,^ ^ ^ ^^(pi+P2-gi-g2)ReA<a:((gi + g2)^) 

2 J 4£/i_C/2 (81) 

Re A^++_ ((gi + ga)', {P2 + QsY, (pi + QsY) 
This expression is directly proportional to ^±(33), so that it will contribute to the imagi- 
nary part of -Boo+- In the next section we will find that the real parts of D]^_^__, D"^^ are 
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7r" 



Figure 3: The 7r"^7r contribution: kinematics in the vr tt CM. 

of the second order in the scattering lengths, so that the contributions of these terms to 
ImA^oo+ are C){a^) and can be neglected. It is convenient to compute the result in the 
CM. of the 7r°7r° pair. With reference to fig. 01 and using eqs. (jSUj) . (|^. we then find 

(Im5oo+),r+7r- =2a^{s3)R+ , 

where 

R^=lj dcos{9) R+ {ip^ + qsf) (82) 

To a good accuracy the integrand is linear in cos (6') — see (JZHJ — so that the average is 
simply the value at ^ = 7r/2. In this approximation we can write 

(ImSoo+).+.- = 2a,(s3)i?+((s)3) (83) 
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where 

{sh = (84) 

Eq. (j83|) reduces to the result in ref. in the hmit S3 4:m'^+. We next consider the 
contribution of intermediate vr^vr^ states (see fig. 12)): 

(ImA^oo+)^o^o = - J ^^^^^^(Pi +P2 - gi - g2)Re A<oo ((^i + ^2)^) ^^^^ 

Re A^oo+ ((P2 + qsf, (Pi + Qsf, {qi + ^2)^) 

We can substantially simplify the computation if we neglect the dependence of Re 
on si and S2. We thus obtain 

(ImA<oo+).o^o = !^!^5^ReA1oo(s3)Re-Moo+(s3) (86) 

o 

We must distinguish the two cases, above and below the tt^tt' threshold. Above the 
threshold we obtain 

(ImAoo+),rO,rO = ^^oo(s3)aoo(s3)i?°(s3), (87) 
(ImSoo+)vr07rO = t;oo(s3)aoo(s3)Re5oo+(s3) (88) 

Below the tt+tt^ threshold the real parts of A^oo+ and A^oo acquire contributions from the 
imaginary parts of -B00+ and of Bqq, see eqs. fTfUjl and (jlZj), and these result in a O^af) 
contribution to Rei?oo+5 

(ReSoo+)vr07rO = -2foo(s3)aa:(s3) [aj.(s3)i?°(s3) + aoo(s3)i?^((s)3)] (89) 

where we have used the results of eqs. flS^ and (jHHj) . The presence of a real part of 5oo+ 
implies an extra contribution to Imi?oo+, eq. (jHHj) . that is however of the third order in 
the scattering lengths, and can be neglected at O(R^). 

Figure m shows the two-pion rescattering in one of the two tt+tt^ channels (the other is 
obtained by the exchange qi ^ 52)- Here the situation is simpler since there is only one 
possible intermediate state: 

(ImA4oo+)7rO,r+ = 2 y 4E E ^^^P'^^P^ - ^2 - qs) ReA^+o ((^2 + qsf) 

Re-Moo+ ((g2 + qsf, {qi+Pz)\ (^1+^2)') ^^^^ 

+ {qi ^ q2) 
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Figure 4: The 7i~^7r^ contribution. 

As anticipated, to a good approximation we can neglect the quadratic terms in eq. (ffij) . 
In this hmit we obtain 



(Im y4oo+)7r07r+ 



(Im4o+)-«-+ + (ImA^ 

+ (si ^ S2) 



.(2) 
^00+ 



(91) 



with 



lom_ I 



(Im A^oV)-»-+ = 2a+o(si)w+o(si)i?" ((s)i - Ai) + ai<o(si)^^ 
where the "velocity" v^q{s) and Aj are defined by 

y+o(s) = 



(92) 



4s,: 



(93) 
(94) 



Finally, we must take into account the effective Svr Stt scattering due to reducible 
diagrams of the type in fig. Ct. By construction, these contributions are at least of Oi^af) 
and at this level of accuracy contribute only to the real part of the amplitude. Following 
the decomposition in eq. (jU^ -ffTO j) . the various rescattering combinations of the type 
in fig. ^ can be divided into two main groups: those which can be reabsorbed into a 
redefinition of Rey4oo+ and those which affect Rei?oo+- We shall start from the latter, 
that are more relevant for the structure of the cusp. 
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The 0{af) corrections to Rei?oo+ arise from diagrams of the type in fig. ^ with the 
identification {na, Hb, vrd = {7r+, 7r°, 7r°} and {tt^', vr^} = {vr^'^Tr^}. We can express all these 
contributions as 

(Re5oo+)fig.[XJ. = -2a^(s3) X [(ImC++_)^+^+ ((5)3) 

+ (ImCj|_)^+^-((s)3)+t;oo((s)3)(ImZ^i!|_)^«^«((^)3)] (95) 

As we shall see in the next section, the three terms in the r.h.s. of eq. correspond 
to the cases where the {vr^,7r^} pair is identified with {7r+,7r+}, {7r+,7r^}, or {7r°,7r°}. 
Summing these three terms we find 

(Re5oo+)fig,[Tl, = -2a^is3)[a+4{s)3)v±i{s)s)R+i{s)3) 

(96) 

where {{s))^ = {Ssq + S3) /A and, given we are already at 0{af), we have neglected the 
tiny P-wave contribution and the difference between Sq and Sq in the {s)i variables. 

Far from the Dalitz plot boundaries, the 0{af) corrections to KeAoo^ could be ignored 
since they can be reabsorbed into a redefinition of R^{s). However, the polynomial form of 
-R°(s) is not appropriate to describe the square- root singularities that occur at the borders 
of the Dalitz plot and, particularly, at tt^tt'^ and tt+tt" thresholds. The latter are described 
at 0(a|) accuracy by the remaining diagrams of the type in fig. The singularities at 
the vr^vr'' threshold, that are obtained with the identification {vr^', tt^} = {ti^tt^}, are 

(SRe Aoo+)ll'^^ " -«oo (-53)^^00 (-53) X (Im A^^Jo)vr+7rO ((-5)3) + (Im aJJJ^+^o ((3)3) 

= -4aoo{s3)voo{s3)a+o{{s)3)v+o{{s)3)R\{{s))3) (97) 

where again we have neglected the tiny P-wave contribution and the difference between 
7r° and vr"^ masses in the {s)i variables. The singularities at the tt+tt^ thresholds, obtained 
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with {tt^', vr^} = {tt+tt"} or {7r''7r+}, are 



{6Re Aoo+)IIj^ ^ -2a+o(si)w+o(si) X Ulm A%(^)^+^o {{s)i) + (ImA+oo) 7,0^0 {{s)i) 

+^±((s)i)(Im 5+00)71+71- + (si ^ S2) 

= -2a+o(si)i;+o(si) [2a+oi{s)i)v^o{{s)i)R%{{s))i) 

+aoo((s)iKo((s)i)i?°((s)i) + 2a,((s)iK_((s)i)i?+((s)i)] 
+ (si^sa) (98) 

In summary, the relevant contribution to the tt^tt^tt^ amphtude (in addition 

to ReAoo^) are: 

Im5oo+ = 2a^{ss)R+i{s);), (99) 
Re5oo+ = -2a^.(s3) [a^(s3)t;oo(s3)i?°(s3) + aoo(s3)t;oo(-53)^^((s)3) 

+a^^i{s)s)v^i{sh)R+i{s),) + 2a+_((s)3)t;±((s)3)i?+(((s))3) 
+a..((s)3)t;oo((s)3)i?°((s)3)] (100) 
Imy4oo+ = aoo(s3)^^oo(s3)-R°(s3) + 2a+o(si)f+o(si)-R° ((s)i - Ai) 

In order to take into account also the 0{af) singularities at the Dalitz-plot boundaries, 
Re A00+ must be modified with the addition of the following extra term: 

Re A00+ ^ R'iss) + <5Re A00+ (102) 

6ReAoo+ = -4aoo(s3)t;oo(s3)a+o((s)3)t;+o((s)3) 

-{2a+o(si)t;+o(si)[2a+o((s)i)t;+o((s)i)i?°(((s))i) 
+aoo((s)i)t;oo((s)i)i?°((s)i) 

+2a,((s)i)t;+_((s)i)i?+((s)i)] + (si ^ S2)} (103) 
4.3 Rescattering in tt+tt+tt^ 

In evaluating the coefficient of the cusp for the vr+vr^vr" decay we need to extract 

from data the real part of the tt+tt+tt^ amplitude. We thus need a suitable param- 

eterization of the latter at the same level of accuracy. Since in the tt+tt+tt^ case the 
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physical region is always above threshold, we do not expect any correction at 0(aj) in the 
decay distribution. This implies that the parametrization (f7H|) for ReC++_ determines 
the real part of the —>■ tt+tt+tt^ decay amplitude at 0(aj) accuracy Since the — >■ 
7r"'"7r"'"7r~ amplitude appears only multiplied by 0(aj) coefficients in the —>■ tt+tt^tt*^ 
rate, knowing the — tt+tt+tt^ amplitude at 0{ai) accuracy is sufficient to the purpose 
of evaluating the cusp effect in Tr+vr'^Tr'^ at the 0(a?) level. 

However, it is worth to stress that the 0(af) corrections to the — > tt+tt+tt^ decay 
amplitude have their own interest: at the border of the Dalitz plot they give rise to square- 
root singularities that could eventually be detected. Their inclusion would therefore 
improve the quality of the tt+tt+tt^ parameterization and could even be used 

to extract an additional information about the vrvr re-scattering. For this reason, in the 
following we shall provide a complete parameterization of the re-scattering effects in the 
—>■ TT+TT+TT^ amplitude to 0{af) accuracy. 

We start from the expression in eq. ()71|) and the parametrization ()75p. where the 
notation of the momenta is defined by 

7r+(gi) + 7r+(g2) + 7r"(g3); Si = (ga + Qsf, etc. 

In analogy with the 71^71^71^ case, we decompose the amplitude isolating explicitly 

the cusp effect related to the vr'^vr'^ ^ tt+tt^ transition. In the physical case, where 
> m^o, this cusp effect is not observable; however, the corresponding amplitude is 
still well defined. As discussed in section this cusp amplitude is more conveniently 
analysed in the unphysical scenario with m^o > m^+. In this scenario the vr^vr'' threshold 
gives rise to two square-root singularities, respectively in si and in S2- To cover the case 
where either si or S2 are below the respective threshold, eq. (ffT|l must be completed as 
follows: 





- = C++- 


- + cK- 


_t.o„(.si) + r>S_ 




si,2 > 4m^o, 


(104) 




- = C++- 




_'t,'„„(s,) + r>S 




si < Amlo, 


(105) 




- = C'++_ 


- + D':i 


(2) 

Voo{si) + iD\[ 


-^'00 (-52) 


S2 < 4m^o 


(106) 



We can choose m^o close to m„+ so that si and S2 cannot simultaneously be below the 
TT^n^ threshold. 
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As far as the two-pion scattering is concerned, we must take into account tt^tt" , 
7T~^TT~ and 7r+7r+ intermediate states. In computing the 7r°7r° contribution to -D^^^_ and 
D^^\__ we will assume that Re A00+ is only a function of S3 as in the parametrization (|7I|l . 
and obtain 

(ImDi_';'l),.vo = a,(si,2)i?°(si,2) (107) 
Analogously, for the tt+tt"'" intermediate state we obtain 

(ImC++_),+^+ =a++{ss)v±{s3)R-^{s3) (108) 

We next consider the contribution of tt+tt^ — > tt+tt^ rescattering, whose general expression 
is. 

(ImA^++_).+.- =i l^^^5\p,+ps-q2~q3)ReM^. {iq2 + QsY) 

ReA^++_ ((g2 + ga)', (gi +P3)', (^1+^2)') 
+ (gi ^ 12) 

Above either of the vr^vr^ thresholds'' in si and S2, we can neglect the contributions of 
the real parts of -D++_ to ReA^++_ that, as we will shortly see, are 0{af) and give a 
0(af) contribution to p09p . Using the parametrization (|7K|l and considering only S-wave 
scattering leads to 

(ImC++_)^+^- = (ImCj|_)^+^- + (ImCS_)^+^- (110) 

with 

(ImCi';'l),+.- = 2a+_(si,2)i^±(si,2)i?+((s)i,2) (HI) 



Note that, although for simplicity of notations we use the same symbol adopted in eq. (j84|l . 

in the — > tt+tt+tt^ case the {s)i variables are defined by 

{sh = (112) 



i.e. we must replace Sq with s'^ with respect to eq. ()84p . If we take into account also the 
tiny P-wave contribution, the above result is modified as follows 

(ImCiVW = 2a+_(si)t;±(si)i?+((^)i) +ai4(si)AV ^4?^ (113) 



^ We recall that we are considering the situation where m.^o > m^i 
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At this point we could consider the (unphysical) case where Si < 4m^o- Here using 

eq. we see that the real part of A4+^ also receives a contribution from the imaginary 

part of -D+l- that, when injected in ()109jl produces a contribution to the real part of 
D^^\__. Another contribution to IieD^^\__ below the tv^tt^ threshold arises from the ImD-i- 
term in Re7W+_, see eqs. (IFk^ and (pij) . Both these contributions are 0(af) and lead to 

iReD^^_U^- =--j -^^5 {p,+ps -q.- qs) ^^^^^ 

= -2v±isi) [R+i{s)i)al{si) + i?°(si)a,(si)a±(si)] (115) 

Finally, we must take into account the effective Svr Svr scattering due to reducible 
diagrams of the type in fig. As in the vt+tt^tt^ case, these contributions are of 

0{aj) and contribute only to the real part of the amplitude. The 0{af) corrections to 
Dt'fL are 



(ReL'5_^|_)gg U], = -a^(si) X ^(Im A^o^+)^o^+((s)i) + (Im A^o^+)^(V+((s)i 

= -4a,(si)a+o((s)i)t;+o((s)i)i?°(((s))i) (116) 

While the remaining 0{af) corrections, that can be absorbed into a redefinition of the 
ReC++_, are 

m^C^+X^l^ = -«++(^3)t^±(^3) X [{lmCXU^-i{s),) + {lmd;l_U^-i{sh) 

+voo{{s)3){lmD''_ll_)^o^o ((5)3) + t;oo((s)3)(ImL'J|_)^o^o ((5)3) 
= -a++{s3)v±{s3) [2a,.((s)3)t;oo((s)3)i?°((s)3) 

+Aa^4{sh)v^{{sh)RH{{s)h)] (117) 

and 

(<5ReC++_)^g^]2^ = -2a+_(si)t;±(si) X [(ImC++_)^+^+((s)i) + (ImCi'|_)^+^- ((s)i) 

+t;oo((s)i)(ImD^'|_),+,-((s)i)] + (si ^ S2) 
= -2a+_(si)t;+_(si) [2a+4(s)i)t;+_((s)i)i?+((s)i) 

+2a+4(.)i)t;+_((.)i)i?+(((s))0 + a,i{s)^)voo{{s)^)R'{{s)^)] 
+ {si ^ S2) (118) 

24 



(1) ^ „ \ , t-t^ /i(2) 



In summary, the relevant contributions to the — > tt^'k'^t: amphtude (in addition 
to ReC++_) are: 



ImL>V+- = a,(si)i?°(si), (119) 
ReL>Ji_ = -2a,(si) [y±(si)i?+((s)i)a,(si)+v±(si)i?°(si)a±(si) 

+2v^o{{s)i)a+o{{s)i)R\{{s)),)\ (120) 
ImC++_ = a++(s3)v±(s3)i?+(s3) + 2a+_(si)v±(si)i?+((s)i) 



(121) 



4.4 The Kl Stt system 



The two Kl — > Stt coupled channels form a system very similar to the one of the two 
— > Stt decays. Similarly to the charged modes, we can decompose the two Kl decay 
amplitudes into regular terms and terms that are singular at the tt+tt" (7r°7r° ) threshold: 

Mqoq = >looo+ J] 5^oW(s^) [e>(si-4m^+)+ie(4m^+ -Si)] , (122) 

i=1...3 

M+_o = C+_o + L'++o^;oo(s3) [e(s3-4m^o) + ie(4m^3 -Si)] ■ (123) 



Concerning the leading amplitudes (ReAooo and ReC_)__o) we shall adopt the following 
phenomenological parametrization: 



Re^ooo(si, S2, S3) = Rl{si, S2, S3) = l + hlj^i^i- OV3<+ 



■8 = 1. ..3 

ReC+_o(si, S2, S3) = i?+(s3) =Al[l+ g+is, - s+J/2ml+ + h+is, - s+J/2mU 



where 



+ 3ml, 



The values of the Kl — > tt+tt^tt" slopes fitted by PDG (that also includes a small term 
proportional to (si -Sa)^) are gl ^ 0.678 ±0.008 and hi ^ /ipdq - {gt/2y = 0.04 ±0.01. 

+ T<^, i Q-rrO ^QOQ Ktiqcv 1 



= m\ + 2ml+ + m^o (/^i 



(Ki ^ SttO) 



(124) 

+ 

(125) 
(126) 



0.008 and ^+ /ipDG - (^+/2)2 = 0.04^ u.u^. 
In the Kl — > Stt^ case the linear term is forbidden by Bose symmetry; the normalization 
of the quadratic term has been chosen such that h\ ^ hpua ~ —(5.0 ± 1.4) x 10~^. 
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Here the visible cusp due to the tt+tt^ tt^ti^ rescattering is expected in the Stt" 
spectrum. The phenomenon is completely analog to what happens in the charged modes; 
however, the relative size of the cusp is smaller because of the inverted hierarchy in the 
leading amplitudes: in the isospin limit A^/A^ = 1/3, to be compared with the analogous 
ratio A'^ /A^ = 2 of the charged modes. 

The calculation of the imaginary parts of the amplitudes (and the real part of the cusp 
coefficient) proceeds exactly as in the charged modes. We report here only the results. In 
the interesting case of the Kl 3tt^ amplitude, the imaginary parts are 

(Im5So),+,- = 2a,{s,)Rl{s,) , (127) 
(Imy4ooo)7r07rO = aoo{si)voo{si)Rl{si,{s)i,{s)i) , (128) 



j=1...3 



with 



{sh = ^^V^ , (129) 
while the 0{af) corrections to the (visible) cusp amplitude are 

(Re5jJo)fig.[lt = -'^voo{si)aa,{si) [a^{si)Rl{si,{s)i,{s)i) + aoo{si)R^{si)] ,(130) 
(Rei?So)fig.[Tl: = -8a,is,)a+oi{s)^Woi{s)^)R-'m)^) . (131) 

For the auxiliary mode, n^n^n^, we find 

(ImL'+_o)7r07rO = a^{s3)Rl{s3, {s) 3,(3)3) , (132) 
(ImC+_o)^±^o = 2a+o{si)v+o{si)Rl ((s)i + Ai^) 

+ a,vUs,)At9l ^ll^ + (^1 - ^2) , (133) 



with 



Ai^ = ^ ^ (134) 



4.5 Decay rates and extraction of the scattering lengths 

The aim of this paper is eminently practical. Our goal is 

• to establish a parametrization of K ^ Sir amplitudes suitable to fit the experimental 
decay distributions at the 10~^ level; 



26 



• describe the cusp effect due to the tt+tt tt^tt" rescattering with a theoretical 
error of a few %. 

Since the cusp effect on the rate is ~ 10%, the two requests are compatible. In this section 
we shall outline the basic strategy for the extraction of the combination of scattering 
lengths a^, as defined in section 1231 from a fit to the tt+tt^tt^ decay distribution. 

From eqs. ()69|) - (f7nj) the differential decay rate for Tr+vr'^Tr'^ is 

|A^oo+r = (ReAoo+ + Re Boo+v±{s3)f + (ImAoo+ + lmBoo+v±{s3)f S3 > 4m^+ 
\Moo+\^ = (Re A00+ - Im Boo+v±{s3)f + (Im A00+ + Re 5oo+f±(s3))^ S3 < 4m^+ 

Expanding the various terms up to 0{af), we can write 

|A<oo+r = (ReAoo+)' + A^ + ?;±(s3)Aeusp + 0(a^) (135) 

where 

Aa = (ImAoo+)' + 4(s3)(Im5oo+)' (136) 
A = / -2ReAoo+Im5oo+ S3 < Aml+ 

1 2ReAoo+Re5oo+ + 2ImAoo+Im5oo+ S3 > ^ml+ ^ ' 

The explicit expressions for the various terms are given in eqs. (|nn|) - ()l()H) . At the same 
level of accuracy, the tt+tt+tt^ decay distribution in the physical region is 

\M++-\^ = (ReC++_)2 + Ac + 0(a3) (138) 

with the 0{af) correction given by 

Ac = 



ImC++_ + t;oo(si)ImDY+- + ^^oo(s2)ImL)^|_ 



+ 2ReC, 



t;oo(si)ReDi.V +t;oo(s2)ReD!,l 



(139) 



and the explicit expressions for the various terms reported in eqs. pi9|l - p2H) . 

The cusp amplitude in eq. p37|) contains both a leading 0(aj) term responsible for 
the negative square- root behavior of the rate below the threshold and an 0{af) term that 
leads to a similar (smaller) behavior also above the threshold (see figure Ej). Both these 
effects are proportional to a^. 
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-0.2 -I 

Figure 5: Illustration of the cusp effects in -k^-k^-k^ . Upper plot: differential decay 

distribution (in arbitrary units, as a function of the vr^vr^ invariant mass) with and without 
the cusp amphtude. Lower plot: relative size of the cusp amplitude with respect to the 
regular term. All plots have been obtained using the values of the scattering lengths and 
their effective ranges from Ref. |^ (see sect. 12 .Sj) . 
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The precision with which the couphng can be extracted from data depends on the 
accuracy of our parametrization of the amphtudes and, in particular, on the theoretical 
expression for Acusp- Since we have neglected 0(af) terms, a priori we should expect 
relative corrections of 0(af) to the value of a^. Given the expected value of the scattering 
lengths and the effect of the 0{af) terms in figure a natural estimate of this error is 
about 5%. A posteriori checks about the size of this error can be obtained by studying 
the stability of the central value of obtained by means of different fitting procedures. 
In particular, it would be interesting to compare the results obtained under the following 
assumptions: 

1. All the a, are treated as free parameters. 

2. All the ai but ax are fixed to their standard values and only is treated as a free 
parameter. 

3. The fit is extended up to the border of the Dalitz Plot with the inclusion of the 
SReAoQ+ term in (frn3|l . 

4. The expressions of R'^'^{s) are modified with the inclusion of cubic terms in (ss — sq) 
and/or quadratic terms in (si — S2). 

5. One of the two Aa,c terms (or both) is ignored [in this way the regular amplitudes 
are re-defined by corrections of O(a^); this, in turn, implies an 0{af) effect on the 
extraction of a^, of the same order of the terms which have not been computed]. 

Finally, it would certainly be quite useful to compare the value of extracted from 
— >■ Svr decays vs. the value extracted in a similar way from Kl Svr decays. 

5 Conclusions 

We have outlined a method that allows to systematically evaluate rescattering effects in 
K Stc decays by means of an expansion in powers of the vrvr scattering lengths. This 
approach is less ambitious than the ordinary loop expansion performed in effective field 
theories, such as CHPT: the scope is not a dynamical calculation of the entire decay 
amplitudes, but a systematical evaluation of the singular terms due to rescattering effects 
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only. In particular, our main goal has been a systematical description of the cusp effect 
in — s> 7r+7r*^7r° P in terms of the tttt scattering lengths. From this point of view, the 
approach we have proposed is more efficient and substantially simpler than the ordinary 
perturbative expansion of CHPT. 

Using this method we have explicitly computed all the 0(af) corrections to the leading 
cusp effect in tx^tt^ti^ , extending the results of Ref. As shown in figure EJ these 

extra terms produce a small square-root behavior also above the tx^tx" singularity. 

The present work allows to reduce the theoretical error on the extraction of oq — 02 
from an experimental analysis of the -k^-k'^ti^ spectrum to about 5%. A similar 

level of theoretical accuracy is also achieved in the case of the —>■ Svr'^ spectrum. This 
level of precision is probably not sufficient to fully exploit the potentially very accurate 
data of NA48, and is also quite above the error on the predictions of ao — 02 in Ref. jlj. 
To reach this level of precision, a complete evaluation of the 0{a^) corrections and — at 
the same time — of the effects due to radiative corrections is needed. 
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